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The semi-inclusive deep inelastic scattering (SIDIS) process is considered. A theoret- 
ical procedure is proposed allowing the direct extraction from the SIDIS data of the first 
moments of the polarized valence distributions and of the first moment difference of the 
light sea quark polarized distributions in the next to leading QCD order. The validity of 
the procedure is confirmed by the respective simulations 
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f 1 dx \d(x) - u(x)] = 0.147 ± 0.039 ^ 0. 
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The extraction of the polarized quark and gluon densities is one of the main tasks of the semi- 
inclusive deep inelastic scattering (SIDIS) experiments with the polarized beam and target. Of 
a special importance for the modern SIDIS experiments are the questions of strange quark and 
gluon contributions to the nucleon spin, and, also the sea quark share as well as the possibility 
Q_i| of broken sea scenario. Indeed, it is known pQ that the unpolarized sea of light quarks is 
o ■ asymmetric, so that the first moments of the unpolarized u and d quark densities do not equal 
52 ■ to each other: 

> 

X 

oj ' The question arises: does the analogous situation occurs in the polarized case, i.e. whether the 
polarized density Am and its first moment 4 Aim = Jq 1 dxAu are respectively equal to Ad and 
Aid or not. 

In ref. [2] the possibility of broken sea scenario was analyzed, considering the results of 
SIDIS experiments on Ag with respect to their consistence with the Bjorken sum rule (BSR) 
predictions. It was shown [2] that using the results of 3\ on the valence quark distributions 
Aiqv obtained in the leading (LO) QCD order, one can immediately estimate the first moment 
difference of the u and d sea quark polarized distributions: 

Ai«-A 1 J= 0.235 ±0.097. (1) 
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At the same time, it was stressed in [2] that this is just a speculation, and, to get the reliable 
results on Aq from the data obtained at the relatively small average Q 2 = 2.5 GeV 2 [3], one 
should apply NLO QCD analysis. The main goal of this paper is to present such a NLO QCD 
procedure allowing the direct extraction of the quantity Aim — Aid from the SIDIS data. 

It is known that the description of semi-inclusive DIS processes turns out to be much more 
complicated in comparison with the inclusive polarized DIS. First, the fragmentation functions 
are involved, for which the information is limited 5 . Second, the extraction of the quark densities 
in NLO QCD order turns out to be rather difficult, since the double convolution products are 
involved. So, to achieve a reliable description of the SIDIS data it is very desirable, on the one 
hand, to exclude from consideration the fragmentation functions, whenever possible, and, on 
the other hand, to try to simplify the NLO considerations as much as possible. 

It is well known (see, for example, [5] and references therein) that within LO QCD approx- 
imation one can completely exclude the fragmentation functions from the expressions for the 
valence quark polarized distributions Aqy through experimentally measured asymmetries. To 
this end, instead of the usual virtual photon asymmetry A^ N = A± N (which is expressed in terms 
of the directly measured asymmetry A\ = (n^ — n^)/(n^+n^) as A\ N = (PbPt j -D) -1 A^ xp ) , 

one has to measure so called "difference asymmetry" A h ^ h which is expressed in terms of the 
respective counting rates 6 as 

Ah-h(„ n 2. x _ 1 (W|| -w^) - (n^ T -nf T ) 
N [ ' P B PTfD{n h u -n\) + {n h n -n^y W 

where the event densities ^jjm) = dN^^^/dz, i.e. n^^^dz are the numbers of events for 
anti-parallel (parallel) orientations of incoming lepton and target nucleon spins for the hadrons 
of type h registered in the interval dz. Quantities Pb and Pt, f and D are the beam and target 
polarizations, dilution and depolarization factors, respectively (for details on these quantities 
see, for example, E] and references therein). Then, the LO theoretical expressions for the 
difference asymmetries look like (see, for example, COMPASS proposal [7j, appendix A) 

+ _ n - AAu v -Ad v _ Au v + Ad v _ 

P 4u V -d V ' d U V + dy ' {) 

A n+-u~ _ 4Ad V ~ Au V A K+-K~ _ Au V A K+-K- _ 
Ady —Uy y Uy 

i.e., on the one hand, they contain only valence quark polarized densities, and, on the other 
hand, have the remarkable property to be free of any fragmentation functions. 

II. Theoretical basis of the procedure 

Let us start NLO consideration with the known jH |H1 13 UD| theoretical expressions for the 
difference asymmetries 

N/h _ N/h 

A%-\x,Q^z)= 9 ^-_ 9 ^ rh (N = P ,n,d), (4) 



5 For discussion of this subject see, for example 0] and references therein. 

6 As usual, one should realize the quantities entering Eq. J2J) not as the pure event densities but 

as the event densities multiplied by the respective luminosities which, in general, do not cancel out - see the 
Appendix. 
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where the semi-inclusive analogs of the structure functions g± and F^, functions g^ h and F^ h , 
are related to the respective polarized and unpolarized semi-inclusive differential cross-sections 
as follows [0] 



d 3 a^n d i(jh Nn _ 4?ra 2 g" /h (x z Q 2 ) (5) 

dxdydz dxdydz Q 2 ' ' ' 



d 3 < 2vra 2 1 + (1 - y) 
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2F 1 N/h (x,z,Q 2 ). (6) 



dxdydz Q 2 y 

The semi-inclusive structure functions g^ n ^ h are given in NLO by 

9,9 Z7r 



+ A^^C^Ee^), (7) 

9,5 



<7i = 9i 



u<-^d 



where the double convolution product is defined as 

[Aq®6C®D)(x,z)= J J^^A q ^5C(x',z')D^y (9) 

v 

The respective expressions for 2Ff <yn ^ h have the form analogous to Eq. (J7J) with the 
substitution Aq — > g, <5C — >■ (7. The expressions for the Wilson coefficients &C qq ( qgm ) and 

Awto^w) = ^(99,9?) + 2 ( X ~ f)^ 1 + ~ y) 2 ) C m{ q9 m) can be found ' for exam P le > in EL Ap- 
pendix C. 

It is remarkable that due to the properties of the fragmentation functions: 

D 1 = Df = Df =Df =Df, 

D 2 = Df = Df =Dl =Df, (10) 

in the differences g^ n — g^ 71 and F^ 71 — Ff^ n (and, therefore, in the asymmetries A 7 ^^ 
and ~ n ) only the contributions containing the Wilson coefficients 5C gq and C qq survive. 
However, even then the system of double integral equations 

o 2 ) ( AAuv - Ady ^ + ^^/( 27r )^®](- p i ~ gg) 

p ,Z) " (4u v -d v )[l + ®a a /(2n)C qq ®](D 1 -D 2 ) ' 

71 (, x >Q 2 'i z ) = (x, Q 2 ; z)\ Uv ^d v 

proposed by E. Christova and E. Leader [I] , is rather difficult to solve directly 7 with respect to 
the local quantities Auy(x, Q 2 ) and Ady(x,Q 2 ). Besides, the range of integration V used in 

7 So that it seems that the only real possibility to extract polarized distributions from the data in NLO QCD 
order is to use some proper fit. At the same time, it is known that a such procedure is rather ambiguous since 
the sea distributions are very sensitive to the choice of initial functions for Aq parametrization (and, especially, 
for Aq parametrization) . 
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0] has a very complicated form, namely: 

x 



— < x < 1 with z < z' < 1, 
x + (1 — x)z 

if x + (1 — x)z >1, and, additionally, range 

x < x' < x/(x + (1 — 

with x(l — x')/{x'{l — x)) < z' < 1 if x+ (1 — x)z < 1. Such enormous complication of the con- 
volution integral range occurs if one introduces (to take into account the target fragmentation 
contributions 8 and to exclude the cross-section singularity problem at Zh = 0) a new hadron 
kinematic variable z = Eh/E]sr(l — x) ('yp cm. frame) instead of the usual semi-inclusive 
variable z h = (Ph)/(Pq) = (labsystem) E h /E^. However, both problems compelling us to 
introduce z, instead of Zh, can be avoided (see, for example 1101) if one, just to neglect the 
target fragmentation, applies a proper kinematical cut Z < Zh < 1, i.e. properly restricts the 
kinematical region covered by the final state hadrons 9 . Then, one can safely use, instead of z, 
the usual variable Zh, which at once makes the integration range T> in the double convolution 
product (9) very simple: x < x' < 1, z^ < z' < 1. Note that in applying the kinematical cut it 
is much more convenient to deal with the total numbers of events (multiplied by the respective 
luminosities - see footnote El and the Appendix) 



iV|\ (TT) (x,Q 2 ) = / dz h rf m) (x,Q 2 ;z h ) (11) 
J z 

within the entire interval Z < z^ < 1 and the respective integral difference asymmetries 10 

a h-h (t0 2, _ i (ivft - jyfi) - (jyft - ivfi) _ 

N {X ^ } z - p B p r f D {N^-N} i ) + {N^-N^)^ [ } 



f l , , N/h N/h^ 
Jz dz h{9l -fli J 

Sldz h {F?' h -F^) 



(N = p,n,d), (13) 



than with the local in Zh quantities n^^^(x,Q 2 ; z h ) and A% h (x, Q 2 ; z h ). So, the expressions 
for the proton and deutron integral difference asymmetries assume the form 11 

^+ 2 _ (AAuy - My) Jj dz h [l + ®f^C w ®](£>i ~ D2) 

P (Au v -d v )tidz h [l + ®l ? ;C qq ®](D l -D 2 ) ' 



A ^- (x Q2) _ (A^y + Ad v ) \l dz h [l + ®%5C qq ®](Di ~ gg) m) 
d {X ^>z- (u v + d v )tidz h [l + ®%C qq ®](D 1 -D 2 ) ' 



8 Then, one should also add the target fragmentation contributions to the right-hand side of Eq. JJJ. 
9 This is just what was done in the HERMES and SMC experiments, where the applied kinematical cut was 
z h > Z = 0.2. 

10 Namely the integral spin symmetries A^ N = dzh g^ N j dzh F[ L N were measured by SMC and HERMES 
experiments (see and also [TUp. 

11 Here one uses the equality g^ h ~ g*y h + 5" which is valid up to corrections of order 0(lod), where 
u>d = 0.05 ± 0.01 is the probability to find deutron in the D-state. 
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where the double convolution product reads 

[Aq ®5C®D] (x, z h ) = I' ^ t % Aq (1) 5C(x', z')D ( ^) . (16) 

J x X J Zh Z \X J \Z J 

With a such simple convolution region, one can apply the well known property of the n-th 
Melin moments M n (f) = Jq 1 dx x n ~ 1 f(x) to split the convolution product into a simple product 
of the Melin moments of the respective functions: 

M n [A ® B] = f 1 dxx n ~ l C —A ( - 1 B(y) = M n (A)M n (B). (17) 
Jo Jx y \y J 

So, applying the first moment to the difference asymmetries A^ + ~ n (x, Q 2 ) and A^ + ~ n (x, Q 2 ] 

given by (JHJ), (fT5j). one gets a system of two equations for AxUy ee Jq 1 da; Auy and Axdy ee 
Jq 1 Ady: 

(4A lUv - A 1 d v )(L 1 - L 2 ) = A e p xp , (18) 



(Ajixy + Aid v )(£i - ^2) = 4 



exp 
d ) 



(19) 



with the solution 



5 L\ — L2 5 L\ — L2 



Here we introduce the notation 

el 



/leap 



AS. 

^d 



exp 



dx a; 



dx Af 



rl OL ~ 

Uu v -d v ) dz h [l + ®—^-C qq ®](D 1 — D 2 ), 

Z Jz Z7T 

/■l a ~ 

(u v + d v ) dz h [l + ^C qq ®]{D l ~ D 2 ), 
z Jz Ztx 



(20) 



(21) 



Li = K 
L 2 = IZ 



T-K 

L d 



Lf 
LI 



TIT 
^d 



(22) 



where 



dz h 



D h q (z h ) 



+ 



2tt 



1 ^ A lC (z')D h q ( Zh 



(23) 



with the coefficient A 1 C(z) = Jq dx 5C qq (x,z). 

Now one may do the last step to get the NLO QCD equation for the extraction of the 
quantity Aim — Aid we are interesting in. Namely, on can use the equivalent of BSR (see |2j 
and references therein for details) rewritten in terms of the valence and sea distributions: 



Axii 



A 1 d= - 
1 2 



9a 
9v 



- {A x u v - Aid v ) 



(24) 
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Using Eqs. (J2~UH2~^) one gets a simple expression for the quantity Aiu— Aid = Jq 1 dx (Au(x, Q 2 ) 
Ad(x,Q 2 )) in terms of experimentally measured quantities, that is valid in NLO QCD : 



Aiu - Aid = - 
2 



9a 



9v 



2A e p xp - 3A e d xp 

nu-L,) ■ (25) 



It is easy to see that all the quantities present in the right-hand side of (25), with the exception 
of the two difference asymmetries Ap + ~ w and A 1 ^^ ^(entering A exp and A d xp , respectively) 

can be extracted from unpolarized 12 semi - inclusive data and can, thus, be considered here 
as a known input. So, the only quantities that have to be measured in polarized semi-inclusive 
DIS are the difference asymmetries A 7 ^ + ~' !T and A\ _7r ^ which, in turn, are just simple 
combinations of the directly measured counting rates. 
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III. Errors on the difference asymmetries 

To check the validity of the proposed procedure let us perform the respective simulations. To 
this end one can use the polarized event generator PEPSI [T^j. First of all, let us clarify the 
very important issue of the errors on the difference asymmetries. At first sight it could seem 
that the difference asymmetries suffer from the much larger errors in comparison with the usual 
asymmetries. Indeed, the approximate formula for the estimation of the statistical error on the 
difference asymmetries reads (see Eqs. (A. 15), (A. 16) in the Appendix) 



^ A il:v - v N : + : N :~ ■ ^ 

So, one can see that, contrary to the usual asymmetries, the difference of the total (for 
both parallel and anti-parallel beam and target polarizations) counting rates for 7r + and 7r~ 
production, N n+ = iV*, + + N++ and N n = NT< + iVj* , occurs in the denominator of the 

expression for 5A^~^ , and it could lead to the large statistical errors on this asymmetry. 
Such situation indeed occurs for the neutron target. However, fortunately, for the proton and 
deutron targets there is an important circumstance which rescues the situation. 

The point is that, unlike the neutron target case (see Fig. 1), the production of positive 
pions on the proton target (see Fig. 2), in the widest region in Bjorken x 13 , essentially exceeds 
the production of negative pions, whereas for the deuteron target the difference in 7r + and 7r~ 
production is not so drastic but is still essential (see Fig. 3). 

It is of importance that though the histograms in Figs. 1-3 are obtained using the 
PEPSI event generator (just summing the events with parallel and anti-parallel beam and 
target polarizations), they represent the general, well established experimentally picture ( 
see section 4.3 and Fig. 10 in ref. [Hj), and peculiar to all known SIDIS event generators 14 . 

12 The spin-independent fragmentation functions D can be taken either from independent measurements of 
e + e _ - annihilation into hadrons or from the hadron production in unpolarized DIS |12| . 

13 This occurs everywhere except for the vicinity of xb = 0, where the ratio N v+ /N 17 approaches unity owing 
to the dominant contribution of the sea quarks |14j . However, let us stress that for the statistical error only the 
difference iV 7r+ -N*~ is of importance, and, at the statistics available to HERMES and COMPASS, iV 7r+ —N n ~ 
is not a small quantity even in the vicinity of the minimal value xb = 0.003 accessible to measurement - see 
below. 

14 Absolutely the same histograms for N^ pd and N£ pd are reproduced using, for example, the unpolarized 
event generator LEPTO |15| . 
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Figure 1: N w and N n obtained with PEPSI for neutron target. The dashed and solid lines 
correspond to 7r + , and n~ productions, respectively. 
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Figure 2: N 77 and N n obtained with PEPSI for proton target. The dashed and solid lines 
correspond to 7r + , and n~ productions, respectively. 



To be sure that concerning the strong asymmetry between N£ d and N£ d the PEPSI event 
generator we deal with strictly follows the real physical picture, one can also compare what 
we have obtained with PEPSI Fig. 4 for the ratio of N£ d and N^ d with Fig. 10 in ref. 



Thus, the differences between the total counting rates N m 



jVfy + JVft and N w ~ 



Ml 



N+< are not small quantities for both proton and deutron targets, and, besides, increase 



TT 



15 Though in the vicinity xb = the ratio N v+ /N 71 approaches unity (see Fig. 4), with the applied statistics 
3-10 6 DIS events (absolutely real for HERMES and COMPASS - see below), the difference N" + -iV 7r ~ entering 
the error significantly differs from zero even near the minimal value xb — 0.003 accessible to measurement. 
Indeed, in the first bin 0.003 < x < 0.006 the PEPSI event generator gives for the proton target N v /N 71 = 
1.085 while - = 2985, and for deutron target N* + /N^ = 1.053 while - = 2020. 
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Figure 3: N 77 and N n obtained with PEPSI for deutron target. The dashed and solid lines 
correspond to 7r + , and n~ productions, respectively. 
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Figure 4: Ratios of N n+ = Nf± + iVf^ and N 71 = iVft + obtained with the polarized 
event generator PEPSI for the different targets. The picture is in good agreement with the 
respective EMC result - Fig. 10 in ref. [Tlj . 

with the statistics. As a result, the respective statistical errors turn out to be quite acceptable. 

Let us illustrate this statement by a simple LO example. Using the GRSV2000LO (broken 
sea) PH] parametrization entering the PEPSI event generator as the input, we generate 3 • 10 6 
DIS events with = 160 GeV 2 (COMPASS kinematics). We then construct the " experimen- 
tal" asymmetries together with their statistical errors using Eqs. (A.l) and (A. 7) from the 
Appendix, respectively. These simulated asymmetries are compared with the theoretical ones 
given by Eqs. (3) - see Fig. 5. One can see that the errors on the simulated asymmetries are 
quite acceptable and that the simulated and theoretical asymmetries are in good agreement 
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Figure 5: Simulated (squares) and theoretical pion difference asymmetries for proton and deu- 
tron targets. The solid lines correspond to the theoretical asymmetries obtained from Eq. (3) 
with GRSV2000LO (broken sea) parametrizations for the valence distributions. 



within the errors. Furthermore, it is seen from Fig. 6 that the extracted valence distributions 
are also in good accordance with the respective input parametrizations. 

IV. Testing of the NLO QCD extraction procedure 



A. Broken sea scenario 



To perform the NLO QCD analysis, we first choose 17 the GRSV2000NLO (broken sea) ^EJ 
parametrization as an input. The conditions of simulations are presented in Table 1 and 
correspond to HERMES and 18 SMC (COMPASS) kinematics. Let us stress that all the cuts 
in Table 1 are the standard cuts applied 19 by SMC, HERMES and COMPASS. The statistics 
3 • 10 6 in Table 1 is the total number of DIS events for both proton and deutron target and for 
both longitudinal polarizations. Since the statistical error on the pion difference asymmetry 

16 Namcly the valence distributions are essential for what follows - see below. 

17 Note that at present the broken sea scenario is argued as the most probable one ^B] (see also discussion 
on this subject at the beginning of this paper and Eq. (I}). 

18 Since the COMPASS muon beam energy 160 GeV is almost the same as that of SMC, 190 GeV, the COM- 
PASS low xb boundary achieved at the asymmetry measurements also should be about the same as the SMC 
0.003. 

19 For example, the important cut on invariant mass W 2 > 10 GeV 2 is applied by these collaborations to 
exclude the events coming from the resonance region. 
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Figure 6: The reconstructed polarized valence distributions (squares). The solid lines corre- 
spond to the respective GRSV2000LO (broken sea) parametrizations. 

depends on N n+ and N n , one need to know the respective semi-inclusive statistics - the total 
(for all available xb) number of pions N^ ot and N£ ot corresponding to 1.5 ■ 10 6 DIS events for 
each target. With the all cuts indicated in Table 1, the PEPSI event generator gives 

^tot \proton — 551281 , Nf 0t \proton = 358654; N^ Qt \deutron = 526747, N^ Qt \deutron — 383826, 

for Ei = 27 GeV, while 

\ pro ton = 582913 , Nl t \ proton = 420709; \ deutron = 559494 , N£ \ deutron = 447599, 

for Ei = 160 GeV. It is of importance that these numbers are absolutely realistic for HERMES 
(in 2000 HERMES already achieved for deutron target iVgJ = 493492, N£ = 402479 - see 
Table 5.4 in ref. [TJj) and even much less than expected by COMPASS (see COMPASS proposal 
0, P- 90 ). 



Table 1: Simulation conditions. A and B correspond to HERMES and SMC (COMPASS) kinematics, respec- 
tively. Here xb and xf are the Bjorken and Feynman x variables, respectively, — is the standard hadronic 
variable and W is the invariant mass of the final hadronic state. 



Kinematics 


Elepton 


x B 


xf 


Zh 


W 2 


Events 


A 


27.5 GeV 


0.023 < x B < 0.6 


x f > 0.1 


z h > Z = 0.2 


W 2 > 10 GeV 2 


3- 10 6 


B 


160 GeV 


0.003 < x B < 0.7 


x F > 0.1 


z h > Z = 0.2 


W 2 > 10 GeV 2 


3- 10 6 



To extract the quantities AiUy, Aidy and, eventually, Aiu — Aid from the simulated 
asymmetries, one should first construct the difference asymmetries together with their statistical 
errors using Eqs. (A.l) and (A. 7) from the Appendix, and then calculate the quantities >ASL 
and Li — L 2 entering Eqs. (20), (25). For A e p xp one should use, instead of integral formula (jzlf 
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the equation (and analogously for A 



exp\ 



a exp 



^bins 

E ^ a p 

i=l 



(xi) (4w y - d v )(xi) / dz h [l + 



Do 



(27) 



where Axi is the i — th bin width. The parametrizations [THj for the fragmentation functions 
and [THj for unpolarized quark distributions are used. Note that here one should not use the 
usual "+" -prescription in the Wilson coefficients C qq , but its generalization, the so-called "A"- 
prescription [201 • The calculation of Li,L 2 is rather simple and can be done using any numerical 
method. 

Let us introduce the additional notation A^qy = f^ max dx Aqy and rewrite 20 BSR in the 
form (24) as 



A x u - A^ 



A{u- A\d 
A{u - A\d 



BSR 



BSR 



9A 



9V 



dx(Auy — Ad\ 



~{A\u v - Ald v ) 



[28) 



(29) 



where x m i n and x max are the boundary points of the available Bjorken x region. It is obvious 
that dealing with the restricted available Bjorken x regions, one can directly extract from 
the measured difference asymmetries namely the quantities Aluy, A\dy and [A\u — A\d}BSRi 
while the "tail" contributions / aw " dxAu v , ^ min dx Ad v and \ $ min dx (Au v - Ad v ) should 
be studied separately, applying the proper extrapolation procedure (see below). 

The results on A^-uy, Aldv and [A*u — AIc[\bsr extracted from the simulated difference 
asymmetries using the presented NLO procedure, are given in Table 2. 



Table 2: GRSV2000NLO (broken sea) parametrization. Results on AJuy, A*dy and [AJu — AIcI]bsr ex- 
tracted from the simulated difference asymmetries applying the proposed NLO procedure. 



Kinematics 


O 2 

mean 


Alu v 


Ald v 


[A$u - Ald] BSR 


A 


2AGeV 2 


0.585 ±0.017 


-0.147 ±0.037 


0.268 ±0.020 


B 


7.0 GeV 2 


0.602 ±0.032 


-0.110 ±0.080 


0.278 ± 0.040 



It is obvious that to be valid, the extraction procedure, being applied to the simulated asym- 
metries should yield results maximally close to the ones obtained directly from the parametriza- 
tion entering the generator as an input. The results for the respective parametrization functions 
integrated over the total < x < 1 region in Bjorken x and over the regions 0.023 < x < 0.6 
(HERMES [3] kinematics) and 0.003 < x < 0.7 (COMPASS kinematics) are presented by the 
Table 3. 

Let us now compare the results from Tables 2 and 3. 

First of all notice that contrary to the actual experiment conditions, the simulations give 
the possibility to check the validity of the extraction method comparing the results of the 
extraction from the simulated asymmetries with an exact answer. Namely, in our case this 
is the integral over the total region of the difference of the parametrizations for Am and Ad 
entering the generator as an input: 

/•0.99 

[Ai« - Aid) exact ~ [Alu - A*d] 2 5 = / dx[Au - Ad] pararnetrization = 0.310, (30) 

J0.0001 

20 As usual, we neglect contributions to Aiq from the unmeasured large xb region 0.6(0.7) < xb < 1 because 
their upper limits given by the unpolarized distributions are very small there - see [13 E]- 
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Table 3: Results on AJ-uy, AJdy, A*u — A{d and [AJu — Ald] B sR obtained from integration of the 
GRSV2000NLO (broken sea) parametrization of the quark distributions over the total and experimentally avail- 
able Bjorken x regions. The fifth column is obtained by direct integration of the respective parametrizations. 
The sixth column is obtained using BSR and the parametrizations for the valence distributions. 



X B 


Q 2 


A\u v 


A\d v 


A\u - A\d 


[A*u - Ald] BSR 


0.0001 < x Bj < 0.99 


2AGeV 2 


0.605 


-0.031 


0.310 


0.315 


0.023 < x Bj < 0.6 


2AGeV 2 


0.569 


-0.114 


0.170 


0.292 


0.0001 < x Bj < 0.99 


7.0 GeV 2 


0.604 


-0.032 


0.309 


0.315 


0.003 < x Bj < 0.7 


7.0 Ge\/ 2 


0.598 


-0.065 


0.262 


0.302 



where symbol [...] nm denotes the n-th line and m-th column of Table 3. 

The second point is that the integral taken directly (without using BSR) from the Am and 
Ad parametrization difference over the available to HERMES region is almost two times less 
than the exact answer (130(1 : 

r0.6 

[Alu - Ald] 35 = / dx[Au - Ad] parametrization = 0.170. (31) 

J0. 023 

This is a direct indication that the HERMES interval in Bjorken x is too narrow 21 to extract 
the quantity Aiu — Aid we are interested in directly. 

However, there is a possibility to avoid this trouble and essentially improve the analysis on 
Aiu — Aid even with the narrow HERMES x B region, applying BSR for Aiu — Aid extraction. 
Indeed, applying Eq. (29) to the HERMES x B region, considering that |<7a/<7v| = 1.2670 ± 
0.0035 and calculating the integrals of the valence quark parametrizations over the region 
0.023 < x < 0.6, one gets 

[A*w - A*rf]f 6 5fl = 0.292, (32) 

and this result (contrary to Eq. (|3ip) is in good agreement with the exact one, Eq. ([30)1 . 

The reason of this good agreement of Eq. (}32|) with the exact answer Eq. dHOJ) is that, 
contrary to the sea distributions, the valence distributions gather far from the low boundary 
x B = (see, for example ^3] and references therein). 

Thus, this exercise with the integrals of the parametrization functions shows that, at least 
within the broken sea scenario, the application of Eq. (|2*4*j) for Aiu — Aid extraction could give 
a reliable result on this quantity even with the narrow HERMES x B region. Namely, one should 
first extract in the accessible x B region the truncated moments of the valence distributions, and 
only then get the quantity [A*w — A*d\BSR applying Eq. (J2HJ)- 

One can also compare elements 55 and 56 from Table 3 corresponding to the SMC (COM- 
PASS) Bjorken x region with the exact answer, element 45 from the Table 3. It is seen that 
even though the integral over the experimentally available region taken directly of the sea 
parametrization difference is now much closer to the exact answer, the application of BSR 
instead of direct extraction significantly improves the situation even for this much wider x B 
region. 

Returning now to the proposed NLO extraction procedure, let us recall that the application 
of BSR in the form (}2*4*j) (see derivation of Eq. (|25|)) is one of the essential elements of the 

21 Note that the proposed NLO extraction procedure has nothing to do with that problem - we just compare 
the integrals of the parametrization over the different Bjorken x regions. 



12 



procedure. Comparing the result of Table 2 on [A*u — A*d]ssR obtained from the simulated 
asymmetries with the HERMES kinematics (element 25 from the Table 2) with both Eqs. (|30p 
and (|32|). one can see that they are in good agreement with each other. Besides, comparing 
the results of Tables 2 and 3 for the COMPASS kinematics, one can see that the results on 
reconstructed A\uy, A\dy and A*u — Aid are in still good agreement with the respective 
quantities obtained by direct integration of the input parametrization over both the total < 
x B < 1 and experimentally available 0.003 < xb < 0.7 regions. 

It is also of importance that even without BSR application, the moments of the valence 
distributions (interesting in themselves) extracted in NLO in the accessible xb regions, are in a 
good agreement with the input parametrization for both HERMES and COMPASS kinematics. 




Figure 7: Polarized valence distributions given by the different parametrizations with the sym- 
metric and weakly broken sea. Solid lines corresponds to the purely symmetric parametrizations 
from refs. [E] and [22], while the dashed line corresponds to parametrization FS2000 (set i- 
in ref. [10J ) with weakly broken sea. All of these parametrizations demonstrate quite similar 
behavior. 



B. Symmetric sea scenario 

Until now we dealt with the broken sea scenario (which seems to be the most probable one - 
see footnote El ) • However, one, certainly, should also investigate an alternative opportunity - 
the symmetric sea scenario. Notice that all the known parametrizations with the symmetric or 
weakly broken 22 polarized sea essentially differ from the only presently known parametrization 

22 The only such parametrization which we know is FS2000 parametrization ^Hj- Ai" — Aid < 0.1 within 
this parametrization, and with respect to the valence distribution it behaves quite analogously to the purely 
symmetric parametrizations - See Fig. 7. At the same time the behavior of the parametrization with the 
strongly broken sea, GRSV2000 (broken sea) , is absolutely different - contrary to the parametrizations with the 
symmetric or weakly broken sea, Ady changes the sign at small xb where the contribution of the sea quarks 
become dominant. 



13 



with the strongly broken sea, GRSV2000 (broken sea), we dealt with previously. However, 
they rather little differ from each other - see Fig. 7. So, for self-consistence, we again choose 
GRSV2000NLO (but with the symmetric sea) parametrization as an alternative input. The 
respective analysis is presented in Table 4. 

Table 4: The upper part presents the results on A*iiy, A\dv and [A^u — AIcI]bsr obtained from integration 
of the GRSV2000NLO (symmetric sea) ^H] parametrization. The lower part presents the results on A*ity, 
A\dv and [A*u — A\d]BSR extracted from the simulated difference asymmetries applying the proposed NLO 
procedure with parametrization GRSV2000NLO (symmetric sea) entering the generator as the input. 





Q 2 


/A\u v 


Ald v 


[Alu - Ald] BSR 


0.023 < x Bj < 0.6 


2AGeV 2 


0.749 


-0.276 


0.121 


0.003 < x Bj < 0.7 


7.0 GeV 2 


0.866 


-0.320 


0.041 


0.0001 < x Bj < 0.99 


2AGeV 2 


0.916 


-0.339 


0.006 


0.0001 < x Bj < 0.99 


7.0 GeV 2 


0.914 


-0.339 


0.007 


Kinematics 


O 2 

mean 


A\u v 


Ald v 


A^u - Aid 


A 


2.4 


0.736 ±0.017 


-0.310 ±0.037 


0.111 ±0.020 


B 


7.0 GeV 2 


0.842 ±0.032 


-0.300 ±0.069 


0.063 ±0.038 



Let us analyze the results from Table 4. First, one can see that for both A and B kinemat- 
ics, the results of reconstruction in the accessible xb region of the all presented in this table 
quantities are again in good agreement with the input parametrization. Thus, the analysis 
performed within the symmetric sea scenario again confirms that the proposed NLO extraction 
procedure satisfies the main criterion of validity - to reconstruct the quark moments in the 
experimentally available xb region. 

On the other hand, performing the reconstruction of the entire quantity Aiu — A\d, one, 
certainly, should not roughly put it to Alu — Aid. It is necessary to carefully estimate the 
unmeasured "tail" | J^ min dx(Auv — Adv) entering Eq. (28), especially dealing with so narrow 
Bjorken x region as the HERMES one. It is clearly seen from the Table 4 where the result on 
A\u — Aid is quite close to the exact answer, zero, for the COMPASS xb region, but in the 
case of HERMES kinematics it indicates rather essential deviation from the zero value. 

C. Low xb uncertainties 

Let us stress that the problem of the unmeasured "tail" estimation is the common and long 
staying problem which, however, in any case should be somehow solved if we wish definitely 
answer the question is the sea symmetric or not. Nowadays the state of art is such that 
the polarized SIDIS experiments use the only method of the low xb contribution estimations 
(see, for example, [Zj): the proper fit to the obtained data on Aq is performed with the 
subsequent extrapolation of the fitting function to unavailable low xb- On the other hand, the 
low xb "tails" of all the existing parametrizations on Aq are obtained using a quite analogous 
procedure. Namely, the parametrization on Aq is extracted in the accessible xb region from the 
fit to the measured inclusive asymmetries and/or structure functions and then is extrapolated 
to low Xb- It is also of importance that the degree of the reliability of the low xb estimations 
applied in the existing parametrizations increases due to that all the parametrizations are 
constructed in the strict accordance with the sum rules on a 3 and a 8 nonsinglet combinations. 
Besides, the constructed parametrizations meet the requirement of agreement with the existing 
DIS data on Y\ and |2U - 
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Table 5: Low-z contributions to \{A\Uv — Aidy) for the different NLO parametrizations. 



NLO parametrization 


jS mm dx(Au v - Ad v )/2 


x min = 0.023, Q 2 = 2AGeV 2 


x min = 0.003, Q 2 = 7.0 GeV 2 


GRSV2000 (broken sea) 


-0.035 


-0.016 


GRSV2000 (symmetric sea) 


0.110 


0.033 


FS2000 (i+) 


0.104 


0.036 


FS2000 (i-) 


0.080 


0.031 


LSS2001 


0.098 


0.032 


AAC2000 


0.116 


0.046 


AAC2003 


0.127 


0.055 



So, we propose to perform the respective estimation of the quantity J^ min dx(Auy — Ady) 
using the maximal number of the latest available NLO parametrizations. The results are 
presented in the Table 5, where the parametrizations from the refs. jTU], JE] and are used. 

Looking at the Table 5 one can conclude that for the HERMES xb region 

1 /-0.023 

dx\{Au v - Ad v )\ <0.13, (33) 



2 Jo 

while for the COMPASS xb region the upper boundary is approximately twice as less: 

1 /-0.003 

-/ dx\{Au v - Ad v )\ < 0.06. (34) 

2 jo 

Notice that one can estimate only absolute value of J^ min dx (Auy — Ady), because we do not 
know which scenario (symmetric or not) is realized in nature. For example, the well known 
broken sea parametrization GRSV2000 gives the negative sign for this quantity while all the 
symmetric sea parametrizations give the positive sign. It is also seen that the restrictions (J3*3j) . 
()34)1 based on the results of Table 5 are rather strong. For safety, we deliberately overestimate 
the upper boundaries choosing the largest numbers from the Table 5 instead of performing the 
averaging procedure over all used parametrizations (just as it was done by SMC [6]). Notice 
also that the restriction ()34j) is consistent with the respective estimation made by SMC (see 
Table 5 in ref. [§]), whereas the upper boundary given by Eq. (|33J) is in four times larger than 
the HERMES estimation 23 . 

Thus, to extract the entire quantity Aiu — Aid we propose to include the upper boundaries 
on | J£ min dx(Auv-Adv) given by the inequalities (|3^j) and (J3~4~j) into the respective systematical 
errors, so that the additional low xb contributions into the systematical errors of HERMES 
and COMPASS look as 

^ WX \ HERMES = ±0.13, (35) 
5 lowx = ±0.06. (36) 

lowx COMpASS V I 

Certainly, Eqs. (|33j) should not be considered as some strict estimations. This is just 
an attempt roughly but with all possible precautions to estimate could HERMES and (or) 
COMPASS under their real conditions answer the question is the sea broken or not. 



2:1 



Notice that it is dangerous to use the simple Regge parametrization for extrapolation in the all rather 



large xb region unavailable to HERMES. The obtained in this way estimations 3 J ° 023 dx Auy — 0.03, 



r-0.023 

lo 



dx Ady — —0.03 seem to be rather underestimated (see discussion on this subject in ref. |2|). 
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V. Discussion and conclusion 



With the (rather overestimated) uncertainties given by Eqs. (JHH|) . (jHEJ) it is quite possible 
that HERMES would not see within the total error that Aim— Aid 7^ 0, if this quantity happens 
too small in reality (for example, about 0.2). However, if this quantity will be about 0.3 (as it 
assumed by GRSV2000 broken sea parametrization) and higher, it could be still possible to see 
this quantity even with the HERMES x b region. On the other hand, it is seen from the Tables 
2, 3, 4 and Eq. (f3T)J) that the COMPASS xb region could allow to catch even small difference 
(if any) between Aiu and Aid. 

In any case, analyzing such a tiny quantity as Aiu — Aid it is very desirable to perform 
a combined analysis with both HERMES and COMPASS data. For example, having at one's 
disposal data on the difference asymmetries in the accessible to HERMES xb region, one 
could involve in the analysis on Aim — Aid the respective COMPASS data from the region 
0.003 < xb < 0.023. The point is that a high statistics, especially at low xb is claimed as one 
of the COMPASS advantages |7j as compared with SMC and HERMES experiments. 

Thus, we have tested the proposed NLO QCD extraction procedure performing the simu- 
lations corresponding to both the broken and symmetric sea scenarios. This analysis confirms 
that the procedure meets the main requirement: to reconstruct the quark moments in the ac- 
cessible to measurement xb region. On the other hand, even with the overestimated low xb 
uncertainty (J3SJ), one can conclude that the question is A\U — Aid equal to zero or not could 
be answered even with the HERMES kinematics in the case of strongly asymmetric polarized 
sea. In any case, the situation is much better with the available to COMPASS xb region. 
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Appendix 

Calculating the asymmetries given by Eq. (fT2"|) together with their statistical errors, one should 
have in mind that, contrary to SMC and COMPASS experiments, in the HERMES conditions 
the quantities Ms ' entering Eq. (JT2J) are not the pure counting rates, but the counting rates 
multiplied 24 by the respective luminosities 3J. Thus, in general case, one should use instead of 
Eq. (H2I) the equation 25 : 



» 7r +_ 7r J- 

A P {n,d) ~ D 



(iVjf - jV T Y)L n - (JVg - Nf~)L n 
(iVf/ - JV£-)L TT + (iVf T + - iVff )L n 



(A.l) 



where luminosities -^TTCTJ.) are defined as 



-^TT(TI) = ( n $)n(U)> ( A - 2 ) 



24 The cancellation of the luminosities is possible only with the special target setup, like the SMC and COM- 
PASS ones 00. 

25 Here the beam and target are assumed to be ideal which means that Pb — Pt = f = 1- Namely this 
assumption is adopted in the PEPSI event generator we use for simulations. 
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where n is th area density of the nucleons in the target and $ is the beam flux. Within the 
paper we do not study the specific peculiarities of the different experimental setups and deal 
only with the event generator where the acceptance a is equal to unity, so that 

-^Ti(TT) = ( n< ^)Ti(TT) acr Ti(tT) ~> ( n *)TI(TT)°"TI(TT) ( A - 3 ) 
Thus, Eq. flA.2|) for luminosities is rewritten in the following, suitable for simulations, form 

L TT(U ) = A r n(U)/( cr TT(U))> ( A - 4 ) 

where iV^^) are the numbers of inclusive events and o"tT(Ti) are ^ ne inclusive cross-sections 
automatically calculated by PEPSI (see ref. ^BJ) for given sets of the kinematic conditions. 

Choosing as the variables in Eq. (|A.1|) the set , N?\ , iV|+ + , and using the general 
formula (see, for example, for the statistical error on the function F of variables xi, x 2 , . . . 

5 2 (F(x u x 2 , ...))= (^-) 5 2 ( Xl ) + (^j 5 2 (x 2 ) + 2 (^j cov( Xl , x 2 ) + . . . (A.5) 

one gets 

6 2 (A;IJ } -) = ^-{(Y- X) 2 L? t [5 2 (<) + 5 2 (N^) - 2cov(<,iV-)] 

+(Y + X) 2 L 2 n [(6 2 N^ + ) + 5 2 (N^~)) - 2 cov(/Vf T + , TVff )} , (A.6) 

where X and Y are the numerator and denominator in the square brackets in Eq. (|A.1J) 
. If the distributions of hadrons N^,N^ ,N^,N^ are Poissonian (low multiplicities n + , 

n~ - see |2Sj) : ^(^n(tl) ^ = \J ^Un) ^ then one can neglect in Eq. (jA.6|) the covariations 
cov(Nfl ,N^~) and cov(iVft + , Nf[) with a result: 

S 2 (A «+-«- ) = ^Mk Wt + - N n) 2 l N n + N n] + « - jVggg + jffl rA 71 

1 P(M) J " ^ 2 [(NftLn + Ntf^-iNftLn + NfrLn)]* ' 1 ' J 

Operating absolutely analogously, one gets for the error on the usual spin asymmetry the 
equation 

MA« + ) - 1 4L T^n(<+<)<< rA R) 

6 1 ^ ] ~ D 2 (Nf+L u + NtfLn)* ' 1 ' 

Notice that namely the equation (|A.8|) for the statistical error on the usual semi-inclusive 
asymmetry was used by HERMES [T7j . 

For the COMPASS experiment = and Eq. (jA.7|) reduces to 

t2( A-K+-ir-\ _ 4 « - iVff ) 2 [< + Nf[] + « - jVf;) 2 [iVf T + + iVff ] 



Let us now choose as the variables the set X\ = iVJ — NJ^ , X 2 = — , so that the 
equation (jA.ljl for asymmetry is rewritten as 



_■ _ 1 XiLtt — XoLti 
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Then Eq. (jA.5|) gives: 

S 2 

4Lf T Lf i X|5 2 (X 1 )+X 1 2 5 2 (X 2 ) 



(A.ll) 



D 2 + A^L-tjJ 

with 

5 2 (X x ) = 6 2 (N^ ) + 5 2 (iVff ) - 2 cov( , JVf ; " ) , (A. 12) 

ff T + ) + 5 2 (iVJf ) - 2 cov(Xf T + , Nf r 



5 2 (x 2 ) = 5 2 (iv T y ) + 5 2 (n:+ ) - 2 cov(iv T y , iVf T ). (A.13) 



Again, with a standard assumption that the distributions of hadrons ■* are Poissonian: 

<5(AyJ^ ^) = \JnJ^I) \ one can neglect in Eqs. (|A.12|) . ()A.13j) the covariations cov(N^ , ) 

and cov(N+, , N*, ). Then Eq. (jA.llj) exactly transforms to Eq. (jA.7j) . 
Let us involve an additional approximation (see [SB], P-7) 

Xi~ X2-Y/2, (A. 14) 

where the quantity 26 Y is defined as Y = (iVfj 4 + Xf T + ) - (JVf^ + iVff) = - A^" . Then 
Eq. (fOT^I reads: 

2/ +— — \ ^ "^^"^TI"^TT 2/~\ ^^^TI^TT 2/ 4- — \ 

5 {A; ^ ] = D~ 2 Y 2 (L n + L n / {Y) = D 2 (L U + L TT ) 4 (N*+ - N*~)* S {N " ^ ~ } 



" (A^+-A^-)2 £2(^ + ^)4- 
With the SMC (COMPASS) target setup L n = L TT and 



(A.15) 



5 2 (A* + -*-) - 1 ^ +N ' (A 1Q) 



It is instructive to reproduce Eq. (|A.16Jl choosing another variables in Eq. (Al): AiV + = 
Nfl - Xf T + , AN 77 ' = Nf[ - Nf, N w+ = Nfl + Xf T + , N*~ = Nf[ + JVff, so that with 
L-^i = Lff (COMPASS target setup) Eq. (|A.1|) for asymmetry reads 

AW+ r - 1 AiV T+ - AN*~ IX . . 

/I — = (A 17 

p(n,d) £) jy7r+ _ ]\fw- D Y' ' 

Then the respective statistical error look as 



2 



s2 ( A ;t:l) = ^f~ 2 Y 2 (AN^ + ) + 5 2 (AN*-) + f^ 2 ^) + 5 2 (iV-))| (A.18) 
' ;i + S(X 7r+ +iV 7r "), (A.19) 



D 2 Y 2 Y 2 



2d 



It is easy to see that Y up to luminosities coincides with denominator Y in the square brackets in Eq. (|A.1|) . 
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where it is again adopted that distributions of hadrons N^,!^. ^ are Poissonian, so that 5(AA^ 



X 




x l 


-x 2 


Y 






+ x 2 



S(N 7T ) = sjNff + Nff, 5(AN* ) = 5(A^" ) = ^Nf[ + JVff and one can neglect cov(AiV 7r , N" 
and cov(AA^~, A^ + ). By virtue of Eq. (TP!) . 

<< 1, (A.20) 

so that one can neglect 27 (X /Y) 2 in Eq. (jA.19)) . Thus, one again arrives at the approximate 
formula (|A.16)1 for the error on the difference asymmetries. 
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